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Abstract 

In this paper, we consider an open tube of diameter e > 0, on the side of which 
a small hole of size e 2 is pierced. The resonances of this tube correspond to the 
eigenvalues of the Laplacian operator with homogeneous Neumann condition on the 
inner surface of the tube and Dirichlet one the open parts of the tube. We show 
that this spectrum converges when e goes to to the spectrum of an explicit one- 
dimensional operator. At a first order of approximation, the limit spectrum describes 
the note produced by a flute, for which one of its holes is open. 

Key WORDS: thin domains, convergence of operators, resonance, mathematics for 
music and acoustic. 

AMS subject classification: 35P15, 35Q99. 

1 Introduction and main result 

In this paper, we obtain a one- dimensional model for the resonances of a tube with a small 
hole pierced on its side. Our arguments are based on recent thin domain techniques of 
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Figure 1: Three different kind of tubes without holes and their resonances. 

[T8] . We show that this kind of techniques applies to the mathematical modelling of music 
instruments. 

Basic facts on wind instruments. 

The acoustic of flutes is a large subject of research for acousticians. Basically, a flute is 
the combination of an exciter which creates a periodic motion (a fipple, a reed etc.) and 
a tube, whose first mode of resonance selects the note produced. Studying the acoustic of 
a flute combines a lot of problems as the influence of the shape of the tube, the study of 
the creation of oscillations by blowing in the fipple. . . see [22], [H], [2], ES] and [5] for nice 
introductions. In this paper, we will not consider the creation of the periodic excitation, 
we rather want to study mathematically the resonances of the tube of the flute and how 
an open hole affects it. Therefore, we simplify the problem by making the following usual 
assumptions: 

- the pressure of the air in the tube follows the wave equation and therefore the resonances 
of the tube are the squareroots of the eigenvalues of the corresponding Laplacian operator. 

- on the inner surface of the tube, the pressure satisfies homogeneous Neumann boundary 
condition. 

- where the tube is open to the exterior, we assume that the pressure is equal to the exterior 
pressure which may be assumed to be zero without loss of generality. 

We can roughly classify the tube of the wind instruments in three different categories, 
depending on which end of the tube is open. See Figure [TJ It is known since a long time 
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that the resonances of the tubes of Figure [T] can be approximated by the spectrum of the 
one-dimensional Laplacian operator on (0, L) with either Dirichlet or Neumann boundary 
conditions, depending on whether the corresponding end is open or not (see for example 
[8]). Notice that this rough approximation can explain simple facts: a tube with a closed 
end sounds an octave lower than an open tube of the same length (enabling for example to 
make shorter organ pipes for low notes) and moreover it produces only the odd harmonics 
(explaining the particular sounds of reed instruments). 

In this article, we study how the one-dimensional limit is affected by opening one of the 
holes of the flute, say a hole at position a e (0, L). At first sight, one may think that it is 
equivalent to cutting the tube at the place of the open hole. In other words, the note is the 
same as the one produced by a tube of length a. This is roughly true for flutes with large 
holes as the modern transverse flute, except that one must add a small correction and the 
length a of the equivalent tube is slightly larger than a. This length a is called the effective 
length. This kind of approximation seems to be the most used one by acousticians. It 
states that the resonances of the tube with an open hole are: a fundamental frequeneyQ 
ir/a and harmonics kn/a, k > 2. However, the approximation of the resonances by the 
ones of a tube of length a is too rough for flutes with small holes as the baroque flute or the 
recorder. In particular, the approximation by effective length fails to explain the following 
observations, for which we refer e.g. to [4] and [26] : 

- the effective length depends on the frequency of the waves in the tube. In other words, 
the harmonics are not exact multiples of the fundamental frequency. 

- closing or opening one of the holes placed after the first open hole of the tube changes 
the note of the flute. This enables to obtain some notes by fork fingering, as it is common 
in baroque flute or recorder. We also enhance that some effects of the baroque flute or of 
the recorder are produced by half-holing, that is that by half opening a hole (some flutes 
have even holes consisting in two small close holes to make half-holing easier). In these 
cases, the effective length a is not only related to the position a of the first open hole, 
which makes the method of approximation by effective length less relevant. 

The purpose of this article is to obtain an explicit one-dimensional mathematical model 
for the flute with a open hole, which could be more relevant in the case of small holes than 
the approximation by effective length. The models used by the acousticians are based 
on the notion of impedance. The model introduced here rather uses the framework of 
differential operators. 

The thin domains techniques. 

The fact that the behaviour of thin three dimensional objects as a rope or a plate can be 
approximated by one- or two-dimensional equations has been known since a long time, see 

1 We use in this article the mathematical habit to identify the frequencies to the eigenvalues of the wave 
operator. To obtain the real frequencies corresponding to the sound of the flute, one has to divide them 
by 2tt 
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and jH] for example. In general, a thin domain problem consists in a partial differential 
equation (E s ) defined in a domain Q e of dimension n, which has k dimensions of negligible 
size with respect to the other n — k dimensions. The aim is then to obtain an approximation 
of the problem by an equation (E) defined in a domain Q of dimension n — k. It seems 
that the first modern rigorous studies of such approximations mostly date back to the late 
80 's: [T5], pQ, [2], [T3], [23]... There exists an enormous quantity of papers dealing with 
thin domain problems of many different types. We refer to (2D] for a presentation of the 
subject and some references. 

In this paper, the domain Q £ is the thin tube of the flute and we hope to model the 
behaviour of the internal air pressure by a one-dimensional equation. It is well known that 
the wave equation in a simple tube can be approximated by the one-dimensional wave 
equation. Even the case of a far more complicated domain squeezed along some dimension 
is well understood, see [T5] and the references therein. We will assume in this paper that 
the open parts of the tube yield a Dirichlet boundary condition for the pressure in the 
tube. In fact, we could study the whole system of a thin tube connected to a large room 
and show that, at a first order of approximation, the effect of the connection with a large 
domain is the same as a the one of a Dirichlet boundary condition, see [3J, [2], [2] and the 
other works related to the "dumbbell shape" model. The main difficulty of our problem 
comes from the different scales: the open hole on the side of the tube is of size e 2 , whereas 
the diameter of the tube is of size e. Thin domains involving different order of thickness 
have been studied in [6], [18], [16], [T7] and the related works. The methods used in this 
paper are mainly based on these last articles of J. Casado-Diaz, M. Luna-Laynez and F. 
Murat. 

Notations and main result. 

For e > 0, we consider the domain 

n £ = (0,1) x (s,0) x (-e/2,e/2) . 

We split any x G Q e as x = {x\, x 2 , x 3 ) = (xi, x). Let a e (0, 1) and 5 > 0. We denote by 
A £ the positive Laplacian operator with the following boundary conditions: 

Dirichlet B.C. u = on (0,e) x {0} x (-e/2,e/2) 

U{l}x(-e,0)x(-e/2,e/2) 
U (a - 5e 2 /2, a + 6e 2 /2) x {0} x (-5e 2 /2, 5e 2 /2) 
Neumann B.C. d v u = elsewhere. 

We denote by Hq(Q, £ ) the Sobolev space corresponding to the above Dirichlet boundary 
conditions. The domain Q e is represented in Figure [2J 

In this paper, we show that, when e goes to 0, the spectrum of the operator A £ converges 
to the one of the one- dimensional operator A, defined by 

D(A) — ► L 2 (0,1) 



A : 



u i — > —u" 
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Figure 2: The domain f2 e . The grey parts correspond to Dirichlet boundary conditions, the 
other ones to Neumann boundary conditions. 

where D(A) = {u G # 2 ((0,a) U (a, 1)) n H%(0, 1) | u'(a + ) - u'(a~) = a5u(a)} and where 
a is the positive constant given by 

a= f |VC| 2 , (1.1) 
Jk 

with ( being the auxiliary function introduced in Proposition 13.21 below. 

Notice that both A e and A are positive definite self-adjoint operators and that 

Vu,veD(A £ ), (A £ u\v) L 2 (Ue) = Vu(x)Vv(x)dx , (1.2) 

Jn e 

Vu, v G D{A) , (Au\v) L 2(q^= I u'(x)v'(x)dx + a5u(a)v(a) . (1.3) 

Jo 

Let < < X 2 £ < \^ < . . . be the eigenvalues of A e and let < A 1 < A 2 < A 3 < ... be 
the ones of A. The purpose of this paper is to prove the following result. 

Theorem 1.1. When e goes to 0, the spectrum of A £ converges to the one of A in the 
sense that 

Vk G N* , A^ >• \ k . 

E e — >0 

Theorem 11.11 yields a new model for the flute, which is discussed in Section [2j The 
proof of Theorem II. II consists in showing lower- and upper-semicontinuity of the spectrum, 
which is done is Sections @] and El respectively. We use scaling techniques consisting in 
focusing to the hole at the place (a, 0,0). These techniques follow the ideas of [T3] (see 
also [16] and [E]). The corresponding technical background is introduced in Section [3j 

Acknowledgements: the interest of the author for the mathematical models of flutes 
started with a question of Brigitte Bidegaray and he discovered the work of J. Casado- 
Dfaz, M. Luna-Laynez and F. Murat following a discussion with Eric Dumas. The author 
also thanks the referee for having reviewed this paper so carefully and so quickly. 
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Figure 3: Left, the first eigenf unction of A, i.e. the fundamental mode of resonance of the 
flute with an open hole. Right, the graphic of the function f a and the intersections with 
the line y = a5 giving the frequencies of the flute. The references values are a = 0.7 and 
a5 = 5. 



2 Discussion 



First, let us compute the frequencies of the flute with an open hole, following the model 
yielded by Theorem 11.11 Theorem 11.11 deals with the spectrum of A e , whereas the res- 
onances of the pressure in a flute follow the wave equation df t u = —A e u (remind that 
A e denotes the positive Laplacian operator). Therefore, the relevant eigenvalues are in 

fact the ones of the operator y ^ J wn i cn are ^"v/^F- Theorem 11.11 shows that 

the frequencies yAf are asymptotically equal to the frequencies \x > such that fi 2 is an 
eigenvalue of A. A straightforward computation shows that /i 2 is an eigenvalue of A, with 
corresponding eigenf unction u, if and only if 

{C sm(/j,x) x G (0, a) 

with some C ^ and with \x > solving 



p -fj.sva.fi 

«<5 = — v ■ , ji 77 := faW , (2.1) 

sm(/ia) sm(/i(l — a)) 



see Figure El 



Using the above computations, we can do several remarks about the resonances of the flute 
with a small open hole, as predicted by our model. 

• The eigenfunctions of A corresponds to the expected profile of the pressure in the 
flute with an open hole, see Figure [3] and the ones of [5], [7] and [26J. 
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• The note of the flute corresponds to the fundamental frequency /x = \/A*. To obtain 
a given note, one can adjust both 5 (the size of the hole) and a (the place of the 
hole). This enables to place smartly the different holes to obtain some notes by 
combining the opening of several holes (fork fingering). We can also compute the 
change of frequency produced by only half opening the hole (half-holing). Notice 
that changing the shape of the hole affects the coefficient a. 

• The overtones of the flute correspond to the other frequencies \x = \/A^ with k > 2. 
We can see in Figure |3jthat they are not exactly harmonic, i.e. they are not multiples 
of the fundamental frequency. This explains why the sound of flutes, which have 
only a small hole opened, is uneven and not as pure as the sound produced by 
a simple tube. In other words, our model directly explains the observation that 
the effective length approximation depends on the considered frequency. Moreover, 
when \i increases, the slope of f a becomes steeper due to the factor \i in (12. ip and the 
solutions of (12. ip are closer to fi = kn. This is consistant with the observation that 
high frequencies are less affected by the presence of the hole than low frequencies, see 
[26] or [25]. However, notice that this is only roughly true since for example one can 
see on Figure [3] that the second overtone is almost equal to 3tt, whereas the fourth 
one is less close to 5tt. This comes from the fact that a = 0.7 is almost a node of the 
mode sin(37nc). 

• Of course, when 5 = 0, we recover the equation sin/i = corresponding to the 
eigenvalue of the open tube without hole. When 5 — > +oo, i.e. when the hole is very 
large, we recover the equations sin(fia) = or sin(/x(l — a)) = 0, which correspond 
to two separated tubes of lengths a and 1 — a (in fact the part (a, 1) is not important 
because this is not the part of the tube which is excited by the Apple). When the 
hole is of intermediate size, the fundamental frequency corresponds to a tube of 
intermediate length a G (a, 1), but the overtones are not the same as the ones of the 
tube of length a. 

• The thin domain techniques used here are general and do not depend on the fact 
that the section of the tube fl £ is a square and not a disk. If the surface g(x) of the 
section of the tube is not constant (think at the end of a clarinet), then the operator 
d%. x in the definition of A must be replaced by -^rd x (g(x)d x .), see [13] . Of course, if 
there are several open holes, then other terms of the type a5u(a)v(a) appear in (II. 3p . 

To conclude, we obtain in this article a mathematical model for the flute with a small 
open hole, which consists in a one-dimensional operator different from a simple Laplacian 
operator. It yields simple explanation of some observations as the fact that the overtones 
are not harmonic. 
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Figure 4: The cube K £ , part of the half-space K = {x G IR 3 , x-i < 0}, and the corresponding 
boundaries. When e goes to 0, the cube K £ converges to K , whereas the hole dK(hole) 
remains unchanged. 

3 Focusing on the hole: the rescaled problem 

When e goes to zero, if one rescales the domain Q £ with a ratio l/(5e 2 ) to focus on the hole, 
then one sees the rescaling domain Q £ converging to the half-space X2 < (see Figure H]). 
The purpose of this section is to introduce the technical background to be able study our 
problem in this rescaled frame. For the reader interested in more details about the Poisson 
problem in unbounded domain, we refer to [21]. 

3.1 The space H\K) 

Let K be the half-space {x G M 3 , x 2 < 0}. For any e > 0, we introduce the cube 

Ke = i) x \w °) x (25? i) 

as shown in Figure HI We denote by dK £ (hole) the part of the boundary (—1/2,1/2) x 
{0} x (—1/2,1/2) corresponding to the hole. We denote dK £ (up) the remaining part of 
the upper face. We also denote by dK(hole) and dK{up) the corresponding parts of the 
boundary of the half-space K. See Figure HI 
We introduce the space H 1 ^) defined by 

H\K) = {ve Hl c (K) , Vv G L 2 (K) and v = on OK (hole)} (3.1) 

and we equip it with the scalar product 

= / V<^.W> . (3.2) 
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We also introduce the space Hq(K) which is the completion of 

C%°(K) = {if E C°°(K) I supp(v?) is compact and (p = on OK (hole)} (3.3) 

with respect to the H 1 scalar product defined in (13. 2p . 

Let x C°°(K) be such that x = 1 outside a compact set, x = on dK(hole) and 
d v \ = on OK (up). Following [21], we get the following results. 

Theorem 3.1. The spaces H 1 (K) and Hq(K) equipped with the scalar product (13 .2p are 
Hilbert spaces and 

H 1 (K) = H l (K)®Rx, (3.4) 

this sum being a direct sum of closed subspaces. 

Moreover, a function u G H l (K) belongs to Hq(K) if and only if it belongs to L 6 (K). 
As a consequence, the splitting of u G H l (K) given by (13. 4p is uniquely determined by 
u = u + ux, where 

1 f 

u = lim — — / u(x)dx 
\K e \ Jk £ 

is the average of u, which is well defined. 

Proof : The direct sum (13. 4p is a particular case of Theorem 2.15 of [21] • The equivalence 
between u G Hq(K) and u G L 6 (K) is given by Theorem 2.8 of [21]. Let u = u + cx with 
u G Hq(K) and c G M. Since u G L 6 (K), we have J K \u\ < C\K e \ 5 ^ 6 and thus the average 
of u is well defined and equal to 0. Since the average of x is we U defined and equal to 1, 
the average of u is also well defined and it is equal to c. □ 



3.2 The function ( 

We now introduce the function (, which is used to define the coefficient a in ( II. ip . 
Proposition 3.2. There is a unique weak solution ( of 

AC = on K 

( = on dK(hole) . . 

d u ( = on dK(up) ^ ' ' 

C = i 

in the sense that ( G H l (K), Q = 1 and 

W^eC^(K) , / VCV<^ = . 

JK 
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Proof : Theorem 13.11 shows that C = X + C with £ G Hq(K). Then, Proposition 13.21 is a 
direct application of Lax-Milgram Theorem to the variational equation 

G H l (K) , / VCV0 = - / A X 4 ■ 

J K J K 

See [21] for a discussion on this kind of variational problems. □ 



The function ( yields a different way to write the scalar product in H 1 ^). 

Proposition 3.3. The function ( is the orthogonal projection of x on the orthogonal space 
ofH^(K) in H X (K). 

Thus, for all u and v in H l {K), there exist two unique functions u and v in Hq(K) 
such that u = ii + uQ and v = v + v(. Moreover, 

{u\v} h^(k) = \ ^u(x).Vv(x) dx + ciu.v , 

J K 

where a is defined by (11.11) . 



3.3 Weak K £ — convergence 

As one can see in Figure HJ if (u £ ) is a sequence of functions defined in Q £ , then the rescaled 
functions w £ are only defined in the box K £ and not in the whole space K. Hence, we have 
to introduce a suitable notion of weak convergence. 

Proposition 3.4. Let {w £ ) £> q be a sequence of functions of H 1 ^^ vanishing on dK £ {hole). 
Assume that 

3C > , Ve > , f \Vw £ \ 2 < C . 

Then, there exists a subsequence (w £n ) ne ^, with e n — > 0, which converges weakly to a 
function u> G H X {K) in the sense that 

G H\K) , / Vw £n V(f > / Vw V(f . 

JK en £ "^° JK 

Moreover, the average of wq is given by 

— 1 f 

w = lim — — | / w £ . (3.6) 



En-S-0 \K F I / /,• 

Before starting to prove Proposition 13.41 we recall Poincare-Wirtinger inequality. 
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Lemma 3.5. (Poincare-Wirtinger inequality) 

There exists a constant C > such that, for any e > and any function <p G H l (K e ), 



\K F 



ip(s)ds 



dx <C 



\V<p{x)\ 2 dx 



(3.7) 



Proof : First, let us set e = 1. The classical Poincare inequality (see [10] for example) 
states that 

1 



ip(x) - 



(p(s)ds 



Ki 



dx<C |V^(a;)rda; . 

'Ki 



Thus, the right-hand side controls the if 1 — norm of <p — Tp. Then, the Sobolev inequalities 
shows that (13.71) holds for e = 1. Now, the crucial point is to notice that the constant C 
in ( 13. 7p is independent of the size of the cube K £ since both sides of the inequality behave 
similarly with respect to scaling. □ 



Proof of Proposition : First notice that Hq(K) is separable due to the density of 
C^°— functions. Hence, H l (K) is also separable and by a diagonal extraction argument, we 
can extract a subsequence e n — > such that for all cp G H 1 (K), \ K Vw £n Vv? converges 
to a limit L(ip) with L(ip) < C\\(p\\^i. By Riesz representation Theorem, there exists 
Wo G H 1 ^) such that L(<p) = (wolv 9 )- 

To prove (I3.6p . we follow the arguments of [TS]. We set w £ = f K w e . Let p G N. 
Lemma l3~5l and the fact that J K \Vw £ \ 2 is bounded, show that there exists a constant C 
independent of e such that f K \w e (x) —w e \ & dx < C. Thus, 

Me < - , / \w e {x) - w £ \ e dx < C . (3.8) 
P Jk 1/p 

By Sobolev inequality, we know that w £ is bounded in L Q (Ki/ p ) (remember that w £ vanishes 
on dKi/ p {hole)). Thus w £ is bounded and up to extracting another subsequence, we can 
assume that w £n converges to some limit /3 G R. By a diagonal extraction argument, 
we can also assume that w £n converges to Wq weakly in L 6 (Ki/ p ), for any p G N. As a 
consequence, (13.81) implies that 

/ \w (x) — (3\ 6 dx < limsup / \w £ (x) — w £ \ 6 dx < C . 

Since the previous estimate is uniform with respect to p G N and since Ku p grows to K 
when p goes to +oo, we obtain that w — (3 belongs to L 6 (K) and thus w — (3x G L 6 (K). 
Theorem 13. II shows that w — fix belongs to Hq(K) i.e. (3 = Wq. □ 
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4 Lower-semicontinuity of the spectrum 



This section is devoted to the following result. 
Proposition 4.1. 



V7c G N* , < limsupA^ < A 



£^0 



Proof : Let {u k ) be a sequence of eigenfunctions of A corresponding to the eigenvalues 
X k . Since A is symmetric, we may assume that (u \vP) ip^i) = for k ^ j. The main idea 
of the proof of Proposition 14.11 is to construct an embedding I £ : Hq(0, 1) —> Hq(Q £ ) such 
that the functions I £ u k are almost L 2 — orthogonal and such that 



7 ' \ ' r-> A fe . (4.1) 

J IWI 2 



The definition of the embedding I £ : H^(0, 1) — >■ ifg(f2 £ ) is as follows. 

Far from the hole: we split the functions u k into two parts u k ^ ^ and u k , ^ , we slightly 
rescale them so that they are defined in (e, a — e/2) and (a + e/2, 1) respectively, and we 
embed both parts in L 2 (Q £ ) by setting 

^ W=U 'Mte75 (a;i - £) ) a ' ld = U "{ a+ l -a- e/2 (Xl ~ °- £/2) 

Near the hole: let ( G ^(K) be as in Proposition 13.21 and let ( = ( + x be the splitting 
given by Theorem 13.11 (where we use that £ = 1 by definition). By the definition of 
Hq(K), there exists a sequence of functions (( £ ) G C£°(K) converging to ( in Hq(K). 
Therefore, there exists a sequence Ce — Ce + X e H ^(K) such that C £ = 1 outside 

a compact set and (C e ) converges strongly to ( when e goes to zero. Notice that we may 
assume that ( £ = 1 outside a compact set of the cube K £ defined in Section [31 We set 
( £ (x) = ( £ ((x - (a, 0, 0))/(fe 2 )) and I £ u k = u k {a)( £ in the cube B £ = (a, 0, 0) + Se 2 K £ . 
Summarizing : the whole embedding I £ is described by Figure El 
Calculing the scalar products: by change of variables, we have 



|2 




Since x is a bounded C°° function and since the volume of K £ is of order 1/e 3 , we have 
Ik W 2 = ^(V^ 3 )- Due to Theorem 13.11 ( £ converges to ( in L 6 (K) and thus 
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u k £ = u k {a) 

Figure 5: The embedding u k = I £ u k of u k G -£^o(0, 1) into H^(Q £ ) (lateral view). 

Therefore, we get that j B \( £ \ 2 = 0(e 3 ). Thus, the L 2 — norm of u k = I £ u k is mostly due 
to the L 2 — norms of (p k and ip k and so, for any k and j, 

(u k \ui) LHn£) = eVl^W) + o(e 2 ) . (4.2) 

On the other hand, we have f B \ V( £ \ 2 = 5e 2 f K |VC £ | 2 - Since (( £ ) converges to ( in H 1 ^) 
and due to the definition (II. ip of a, f K |VC £ | 2 converges to a. Therefore, for any k and j, 

Vip k V(fi+ I Vi) k Vi)l+u k {a)u ] {a) [ |VC| 2 



/ Vu k Vui = [ Vy^V^ + [ 



Us 



d x u k (x)d x u j (x)dx + / d x u k (x)d x u j (x)dx 

J a 

+ 6u k (a)u>(a) j |VC £ | 2 J + o(e 2 ) 
d x u k (x)d x u j (x)dx + adu k (a)u j (a) ) + o(e 2 



e 2 (Au k \u j ) L 2 + o(e i 



(4.3) 



Hence the previous estimates yield the limit (14. ip . 

Applying the Min-Max formula: for e small enough, ( 14. 2 p implies that the functions u k are 
linearly independent. Due to the Min-Max Principle (see [21] for example), we know that 



\ k < 



mm max 







2 




| J2i=i cJ £ uP' | 


2 



(4.4) 



The above estimates (14. 2 p and (14 .3p show that, for any c 6 M. k , we have 



+ o(l) 







l 2 _ (AYLw? 1 






IEti c ^ MPi l 


2 II V fc 

II £ji=i 


CiUP*\\ 2 L2 
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where the remainder o(l) is uniform with respect to c when e goes to zero. Using the 
Min-Max Principle another time, we get 

mm max — - r = A + o(l) . 



This finishes the proof of Proposition I4.ll □ 



5 Upper-semicontinuity of the spectrum 

Let e > and let (u k ) be a sequence of eigenfunctions of A £ corresponding to the eigenvalues 
(A*). We can assume that the functions u k are orthogonal in L 2 (Q £ ) and that ||i4Hl 2 = £• 
To work on a fixed domain, we set f2 = (0, l) 3 and we introduce the functions v k = Ju k 
where J is the canonical embedding of H l (Q e ) into that is that 

Ju k {y) = v k (y) = v k ( Vl ,y) = u k ( yi ,ey) . 

We have 

+ AW and Itfll^l. 

By multiplying the previous equation by v k and integrating, we get 

/ \d yi v k £ ? + y 2 \dyv k e \ 2 = \ k e . (5.1) 
Jn e 

Proposition 14.11 shows that (X k ) £ >o is bounded. Therefore, up to extracting a subsequence, 
we may assume that (A*) converges to Aq = lim inf £ _> X k when e goes to and that (i;*) 
converges to a function v k £ H l {VL), strongly in iJ 3//4 (f2) and weakly in H l (VL). Moreover, 
(15. ip shows that v k depends only on y\. In the following, we will abusively denote by v k , 
either the function in if 1 (fi) or the one-dimensional function in i? 1 (0, 1). 

The purpose of this section is to use the methods of [18] (see also [16] and [IT]) to prove 
the following result. 

Proposition 5.1. For all k £ N* ; the function Vq is an eigenfunction of A for the eigen- 
value Aq. 

Proposition 15.11 finishes the proof of Theorem 11.11 since we immediately get the upper- 
semicontinuity of the spectrum. 

Corollary 5.2. 

VA;£N*, liminfA^>A fc . 
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Proof : We recall that the functions v £ are orthonormalised in L 2 (Q) and converge strongly 
in L 2 (Q) to Vq. Thus, the functions are also orthonormalised. Since Aq = liminf e _^ A*, 
we know that Aq < Aq < . . . < Aq. Then, Proposition 15.11 shows that Aq, . . ., Aq are k 
eigenvalues of A with linearly independent eigenfunctions, and thus that the largest one 
Aq is larger than X k . □ 



The proof of Proposition 15.11 splits into several lemmas. To simplify the notations, we 
will omit the exponent k in the remaining part of this section and we will write u £ for u £ , 
v £ for v £ etc. 

Lemma 5.3. Let B £ C £l £ be any cube of size e and let T £ be one of its faces. Then, 

— J u £ {x)dx = — J u £ {x)dx + o(l) . (5.2) 

As a consequence, vq satisfies both Dirichlet boundary conditions t>o(0) = t>o(l) = 0. 

Proof : We split the cube in slices B £ = U se [ 0t£ ]T £ (s) with T £ = T £ (0) and we set x = (s, x) 
with x e Tffs). For each s, we have 



u £ (s,x)dx — / u £ (0,x)dx 
v e {s) Jr s (o) 




Vm £ (^, x) \ d£dx 



|Vm £ (^, x)\ 2 d^dx . 



To show (15 .2p . we integrate the above inequality from s = to s = e and we notice that 
||Vu e || L 2 = A £ ||u e || L 2 = X £ e = 0(e). 

The fact that Vq(1) = follows from v £ (l,x) = and the strong convergence of v £ to Vq 
in if 3 / 4 (f2). To obtain the other Dirichlet boundary condition, we apply (15.21) to the cube 
B £ = [0,s] x [—£,0] x [— e/2,e/2] at the left-end of fl £ . Since u £ vanishes on the upper 
face of B £ , the average of u £ goes to zero in B £ . Applying (15. 2p again, the average of u £ 
goes to zero on the left face of B £ . Thus, the average of v £ goes to zero on the left face 
r = {0} x [—1,0] x [—1/2, 1/2] of Q and hence J r v (0,y)dy = because v £ converges to 
v in i7 3 / 4 (f2). Since v does not depend on y, this yields v (0) =0. □ 



We now focus on what happens close to the hole at (a, 0,0). To this end, we use the 
notations of Section [3] and we introduce the functions w £ e H l (K £ ) defined by 

Vx G K £ , w £ (x) = u £ ((a, 0, 0) + 5e 2 x) . 
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The functions w £ will be useful to study the behaviour of u £ in the cube 

B e = (a- e/2, a + e/2) x (-£, 0) x (-e/2, e/2) . 
We show that they weakly converges to t>o(a)C m H l (K) in the following sense. 
Lemma 5.4. For a// ip G H l (K), 

Vw e Vip > v (a) / VCVy? • 

K e Jk 

Proof : We have 

j_ /■,_.,. i i . y . A 



VwJ 2 = — / |VuJ 2 < / IVuJ 2 = —A, / \u ' 2 7 ; 



Moreover, the average of w £ in iT £ is equal to the one of u £ in B £ , which converges to vo(a) 
due to Lemma [5.31 and the convergence of v £ to vq in if 3//4 (f2). Applying Proposition 13. A\ 
we obtain the weak convergence of a subsequence of w £ to a limit wq, whose average is 
Wq = Vo(a). To prove Lemma [5.41 it remains to show that wq = Vo(a)(, which does not 
depend on the chosen subsequence (e n ). 

Let <f G C£°(K) and assume that e is small enough such that supp((/>) C K £ . We set 



Vx G B £ , (p E {x) = <p 



x — (a, 0, 0) 



8e 



2 



and we extend ip £ by zero in Q £ . Since 

\\u £ \\ L 2=e and ||<p e ||£2(n„) = S 3/2 e 3 \\cp\\ L 2 {K ) , 

we get 

Vw £ Vy? = t"? / Wu £ W(p £ — ~Y~n [ ^u £ (p £ = y^t [ u £ ip £ > 



&2 _ 6 r£ fe2 ™ cir6 fe2 ^ £ 



Thus, ifo is orthogonal to Cq'(K) and hence to Hq{K) and Proposition 13.31 implies that 
— ^oC- Since we already know that Wq = Vq(o), Lemma [5.41 is proved. □ 



Proof of Proposition \5.1\ : We have shown in Lemma 15.31 that Vq satisfies Dirichlet 
boundary condition at X\ = and x\ = 1. Let (ft G Hq(0, 1) be a test function. We also 
denote by (ft the canonical embedding of (ft into We embed (ft into Q £ by setting 

(ft £ = I £ (ft, where I £ is the embedding introduced in the proof of Proposition 14.11 Using the 
notations of Figure El we have 

/ Vu £ V(ft £ = I Vu £ Vp £ + (ft{a) I Vw £ VC £ + / Vm £ V^ £ (5.3) 

Jn e JxxKa-e/2 J B E Jx 1 >a+e/2 
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The limits of the different terms are as follows. First, notice that 

Vn e V0 e = A e / u £ <p £ = e 2 \ £ / v £ J<p £ 
n E Ju e Jn 

where J<p £ is the canonical embedding of £ in H 1 ^). Obviously, J<p £ converges to J<p in 
L 2 {VL) and we know that v £ converges to v in L 2 (Q). Thus, 

Vu £ V<p £ = e 2 X [ v 4> + o(e 2 ) = e 2 X [ v (p + o{e 2 ) . 
n £ Jn Jo 

In the parts x\ < a — e/2 and x% > a + e/2, we know that v £ converges to vq weakly in 
H l (Q) and obviously J(p e and Jip £ converge to <fi strongly in H l . Moreover, notice that 
Jtp £ and J%l) e only depends on x\. Hence, 



Vu £ Vif £ + / Vm £ V^ £ 

xi<a— e/2 Jx±>a+e/2 



e 2 



J2 



( / d xl v £ d Xl {J(f £ ) + / d Xl v £ d Xl {Jip £ )) 

\J x 1 <a-e/2 Jxi>a+e/2 J 

J d Xl v d Xl (p + J d xi v d xi (f) \ +o(e 2 ) . 



= e 2 I d Xl v d Xl 4> + o(e 2 ) . 
Jo 

The term of (15. 3p in the box B £ is more delicate, but all the work has already been done 



in Lemma [5.41 Indeed we have 



Vu £ V( £ = 6e z / Vw £ V( £ . 

B e JK E 

By definition ( £ converges to ( strongly in H l (K). Thus, Lemma [5.41 implies that 



Vw e VCe = Se 2 v (a) / VCVC + o(e 2 ) = aSv (a)e 2 + o(e 2 ) . 

B e J K 

In conclusion, when e goes to 0, Equality (15.31) shows that 

A / v <f)= / d Xl v d xl <p + a5v {a)<f)(a) . 



Since this holds for all <p £ Hq(0, 1), going back to the variational form of A given in (11.31) . 
this shows that Vq is an eigenfunction of A for the eigenvalue Ao (remember that ||fo||,L 2 = 1 
and so vq is not zero). □ 
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